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I  INTRODUCTION 

Photoinduced  Bragg  Grating  in  optical  fibers  were  first  demonstrated  by  Hill  et  al7 
in  1978.  In  1988  Meltz,  et  al.^  developed  an  external  Bragg  grating  writing  technique 
which  had  more  flexibility  in  the  choice  of  the  Bragg  period.  Since  then  and  due  to  their 
unique  properties  a  number  of  sensing  concepts  have  been  proposed  and  demonstrated. 
Most  of  these  sensors  are  configured  as  back  reflectors.  A  more  general  configuration,  a 
Tapped  Bragg  Grating  (TBG)  (see  Fig.  1),  was  first  demonstrated  by  Meltz  et  al.^  in 
1991.  A  complete  theory  for  the  fabrication  and  operation  of  TBG’ s  is  still  incomplete. 

In  this  paper  a  simple  model  of  Tapped  Bragg  Gratings  is  developed.  This  model 
will  allow  the  understanding  of  the  basic  parameters  necessary  to  fabricate  TBG’s. 

Several  questions  have  to  be  answered  to  effectively  fabricate  a  TBG  sensor.  Some  of 
these  questions  are.  What  must  the  Bragg  angle  and  period  (0b  and  Ab)  be  in  order  to 
effectively  scatter  radiation  out  of  the  fiber  at  a  desired  scattering  angle  0  (refer  to  Fig.  1 
for  a  graphical  definition  of  the  parameters)?.  How  is  the  energy  spread  out  of  the  fiber 
for  a  given  0b  and  Ab?  and  what  is  the  band  pass  of  the  scattered  radiation  for  a  given  0b 
and  Ab  ?  (this  information  is  important  to  determine  how  many  sensors  can  be  placed  in  a 
single  fiber). 

In  section  II  the  model  and  the  TBG  fabrication  parameters  will  be  described. 

This  is  the  largest  section  and  will  be  subdivided  into  several  smaller  sections.  This 
subsections  include  “The  Induced  Current  Approximation”,  “Induced  Current  Expansion”, 
“Vector  Potential  Calculation”.  The  band  pass  of  a  Bragg  grating  and  the  band  pass  of  a 
TBG  Avill  be  calculated  for  two  different  fiber  modes  in  section  III.  Section  IV  will 
calculate  the  Poynting  Vector  and  Scattered  Power.  Section  V  will  show  how  the  energy 
is  spread  out  of  the  fiber  for  two  different  fiber  modes.  Finally  in  section  VI  the  summary 
and  conclusions  will  be  given. 
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Fig.  1 .  Basic  parameters  used  in  the  model 


II  TBG  MODEL 

To  answer  the  questions  addressed  in  the  introduction,  the  induced  current 
approximation^  will  be  used  to  model  the  effects  of  a  TBG  in  an  optical  fiber.  This  section 
will  be  devided  into  three  subsections.  Section  Ila  will  introduce  the  induced  current 
approximation.  Section  Ilb  will  expand  the  current  into  its  cartesian  coordinates.  Section 
lie  will  calculate  the  vector  potential  from  where  the  scattering  fields  will  be  determined. 

Ila.  The  Induced  Current  Approximation 

In  the  induced  current  approximation,  radiation  fields  produced  by  perturbations  to 
the  index  of  refraction  of  the  core  of  optical  fibers  are  assumed  to  originate  from  induced 
currents  in  the  fiber.  These  currents  take  the  form  of  ^ 


J  =  i^  k-(n^-ny-E  (1) 

where  no^-n^  is  the  square  of  the  differences  of  the  perturbed  to  unperturbed  index  of 
refraction,  E  is  the  unperturbed  electric  field  inside  the  core  of  the  fiber,  and  k  is  the  wave 
number.  Assuming  that  the  fiber  axis  is  along  the  z-direction  and  that  Kb  is  in  the  x-z 
plane,  then  the  quantity  no^-n^  in  Eq.  1  can  be  modeled  by  a  simple  trigonometric  function 


n^  -  n^  =  2noAn  •  cos(Kg  •  f') 

=  2n„An-  0(R-r)-0(L/2-|z|)cos{xKB  cos(0b)  -  zKBsin(0B)} 


(2) 


where  An  is  the  difference  in  the  index  of  refraction  of  the  perturbed  to  unperturbed 
regions  in  the  fiber  core.  The  step  functions  in  Eq.  2  are  used  to  define  the  extent  of  the 
TBG.  The  argument  of  the  cosine  is  the  dot  product  of  the  Bragg  wave  number 
(Kb=27c/Ab)  times  the  position  vector  in  Cartesian  coordinates.  Figure  1  shows  a 
graphical  representation  of  the  parameters  used  in  this  model.  The  electric  field  in  Eq.  1 
will  depend  on  the  mode  traveling  in  the  fiber.  There  are  two  modes  that  will  be 
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considered  in  this  paper,  hybrid  modes  and  transverse  modes.  For  a  hybrid  field  mode  the 
electric  is  given  by’ 


E(r,(l),z)  =  e(r,<|))  •  e‘P" 


e  =  e,(r,(j))  •  e,  +  e^(r,<|))  •  e^  +  e,(r,(t>)  •  e. 


where 


e,(r,(|))  = 


^2'^  v+l(^  J^) 

Jv(u) 


fv(4>) 


a,Jv-i(u:^)-a2J,„(u:^) 
e,(r,<|))  = - 

Jv(u) 

iu  Jv(U^) 

ez(r,<|))  =  .  fvW 


RP  J,(u) 


fv(4>)  = 

gv  (<!>)  = 


[cos(v(l)),  for  even  v 
[sin(v(|>),  for  oddv 
[-  sin( v(|)),  for  even  v 
I  cos(v(l)),  for  odd  v 


a,  = 


and 


F-1 

2 

F  +  1 


(3) 


(4) 


F  = 


uw/  b,  +b 


•  b  1  fJv-i(u)  Jv^i(u)l.  1  [k,_i(w)  ^  K^,,(w)] 

'  2u[j,(u)  J,(u)  J’  ’  2w|k,(w)  K,(w)  J 


(5) 


u  =  R(k’nL-p’y'’  ;  w  =  R(3’-k’n’y'^  ;  v  =  Rk(nL -n’,)’'’ 


For  a  transverse  field  mode  the  electric  field  has  only  one  component  and  is  given 


by 


E  = 


J.(u-) 

R  e'P^-e. 


J,(u) 


(6) 


The  induced  current  approximation,  therefore,  assumes  that  a  Bragg  grating  is 
formed  by  a  large  coherent  collection  (  N  =  L/Ab)  of  identical  radiators  and  does  not  take 
into  account  multiple  scattering  effects  between  them.  The  first  approximation  (identical 
radiators)  breaks  down  for  long  gratings  since  the  intensity  of  the  electric  field  will  drop  as 
energy  is  scattered  out  of  the  fiber.  The  second  approximation  (no  multiple  scattering 
effects)  breaks  down  as  the  Bragg  angle  approaches  180°  (the  so-called  back  reflection 
mode).  In  the  limit  of  small  reflection  coefficient  for  each  Bragg  line  (which  is  the  case 
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for  most  fiber  Bragg  gratings)  only  the  first  reflection  will  be  significant  and  the  second 
approxdimation  will  remain  valid  . 

Ilb.  Induced  Current  Expansion 

All  of  the  calculations  in  this  effort  will  be  performed  in  a  Cartesian  coordinates 
system.  Therefore,  by  transforming  the  induced  current  (Eq.  1)  into  this  coordinate 
system  it  is  simple  to  show  that 


=  a  •  E^  =  a  •  e,  cos((t))  -  a  •  e^  sin((l)) 

Jy  =  a  •  E,^  =  a  ■  e,  sin((|))  +  a  •  e^  cos((t))  (7) 

j  zz  a  *  e 

z  z 

where,  for  clarity,  all  the  parameters  (except  for  the  electric  field  components)  have  been 
absorbed  in  the  parameter  “a”.  By  separating  the  “even”  and  “odd”  terms,  it  can  be  easily 
shown  that  for  the  hybrid  modes 


Jev'en  _ 

X  =-a 


Jv-.(u-) 

Jv(u) 


cos(l  -  v)(|)  +  aj 


Jv.,(u-) 

Jv(u) 


COS(l  +  v)(j) 


^iPz 


=  -a 


Jv-,(u^) 

-ai  — — — sm(l  -  v)(t)  +  a, 
Jv(u) 


Jv.i(u-) 

Jv(u) 


sin(l  +  v)(j) 


(8) 


Jy'“  =  -a 


J..-,(u-) 

Jv.(u) 


sin(l-  v)^  +  a2 


J..,(u-) 

J.(u) 


sin(l  +  y)^ 


rodd 


a,  -  cos(l  -  v)^  -  a^  -  cos(l  +  y)^ 

Jv(u)  J,(u) 


j  iu  [cos(v<|))]  even 

^  RP  J,(u)  [sin(v(l))J  ®  odd 


(9) 


and  for  the  transverse  modes 
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h 

h 


-  _a - R_  sin  <j)  ■ 

J,(u) 

=  +a - ^coscj)  • 

J,(u) 


(10) 


J,=0 

lie.  Vector  Potential  Calculation 

From  the  analytical  expressions  of  the  induced  currents  (Eqs.  8  - 10),  the  radiation 
fields  can  be  calculated  using  standard  antenna  theory,  that  is,  by  calculating  the  vector 
potential 

fj(r',<|)',z')-e'“‘*''“^’‘'^-dV'  =  — •e““"*-M  (11) 

47ts  47CS 

where  s  indicates  the  scattering  direction  and 

s'cos(x')=  r'sin(0)cos(<j>-<j)')  +  z'-cos(0)  (12) 

where  the  primed  coordinates  refer  to  a  “cylindrical”  coordinate  system  describing  the 
fields  inside  the  core  of  the  fiber,  while  the  unprimed  coordinates  refer  to  a  “spherical” 
coordinate  system  describing  the  scattered  fields  outside  the  fiber.  Fig.  2  shows  the 
relationship  between  all  coordinates. 


Fig.  2:  Coordinate  system  used,  s  indicates 
the  scattering  direction  while  s’  indicates 
a  position  inside  the  fiber. 
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To  calculate  the  vector  potential  A,  the  integral  in  Eq.  11  will  be  calculated 
sequentially  by  first  integrating  over  the  variable  z’,  then  by  integrating  over  the  variable 
(j)’  and  finally  by  integrating  over  r’. 

Integral  over  the  variable  z’ 

The  only  terms  containing  the  variable  z’  are  in  the  index  of  refi'action  (Eq.  2)  and 
in  the  exponential  factor  exp(ipz’).  These  terms  are  common  to  all  the  components  of  the 
vector  J,  therefore  the  intagral  that  must  be  solved  is 

L/2 

I,  =  j cos(r'  cos((t)')  •  Kg  •  cos(0b)  -  z'-Kg  •  sin(0g))  •  e‘P^'  •  e-^'"“<«>^'dz'  (13) 

-L/2 


By  expanding  the  first  cosine  in  terms  of  exponentials  and  by  introducing  a  new  set  of 
dummy  parameters  a,  b,  and  d  given  by  r’cos((t)‘)  Kbcos(0‘)  ->a ,  Kb  sin(0B)  b  and 
kncicos(0)->d)  it  can  be  easily  shown  that 


L  =  r  •  dz'  =  e“ 

2  L 


siiH^(-b  +  p-d)- 
2i(-b  +  p-d) 


■  +  e 


sin|^(b  +  p  -  d) '  - 
^i(b  +  p-d) 


_  ^ir’KBCOs(0B)cos(<|)’)  . 


sid  (~Kg  sin(0g )  +  P  -  kn^j  cos(0))  • 


2J 


-Kg  sin(03 )  -f  p  -  kn^i  cos(0) 


-  + 


sm 


-f  e 


■ir'KB  cos(6b  )cos(4)’) 


(Kg  sin(0g )  +  P  -  kn^,  cos(0))  • 


Kg  sin(0g )  +  P  -  kn^,  cos(0) 

Since  the  bound  modes  are  always  such  that 

k  •  n„  <  p  <  k  •  n,„ 

only  the  first  term  will  give  an  appreciable  contribution.  Therefore,  Ii  reduces  to 


I.  =  e 


sinl 

ir'Kg  cos(0B )  cos((|)’) 


(-Kg  sin(03)  +  p  -  kn^i  cos(0)) 


2J 


-Kg  sin(0g)  -I-  p  -  kn^i  cos(0) 


(14) 


(15) 
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If  the  length  L  of  the  TBG  is  large,  then  Ii  will  be  appreciably  different  from  zero 
only  when  the  argument  of  the  sine  vanishes.  From  this  one  obtains  the  Bragg  condition. 

First  Bragg  Condition  ( 1 6) 

This  is  one  of  two  equations  needed  to  properly  fabricate  a  TBG.  It  can  be  easily  shown 
that  this  condition  is  simply  a  statement  of  momentum  conservation  along  the  z-direction. 

Integral  over  the  variable  (f)’ 

Now  we  will  calculate  the  second  integral,  i.e.,  we  will  integrate  on  the  azimuthal 
variable  <j)’.  Since  each  component  of  the  vector  J  has  a  different  dependence  on  <|)’,  all 
components  will  have  to  be  calculated  separately.  We  will  first  calculate  the  x  component 
for  the  even  hybrid  modes.  The  other  components  and  modes  can  be  calculated  similarly. 


and 
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a  =  r'  Kg  cos(0g ) 
b  =  r'kn^,  sin(0) 


(19) 


in  a  similar  way,  it  can  be  shown  that  for  the  other  components 


M»<“  =  271  •  i'-' 


Jv-l(u-)  Jv.l(u-) 

h  ,  .  f  •  Jv-i(Q)sin(v  -  l)v|/  +  a,  •  J^,,(Q)sin(v  +  l)vi/ 

JvW  Jv(u) 


=  27t  •  i 

y 


v-l 


T  r  r  •  J  v-l  (Q)  sin( V  -  1) V  -  a^  -  ^  ^  •  J ,,,  (Q)  sin( v  +  l)v|/ 

Jv(u)  Jv(u) 


Mf  =  2k  ■  i''-’ 


Jv-l(u-)  Jv.i(u-) 

— j— - —  Jv-,(Q)cos(v- 1)1(/  +  a, . •  J,„(Q)cos(v+l)\|/ 

Jv(u)  Jv(u) 


^  iu  f-i''Jv(Q)-cos(vv|/)  V  =  even 

RP  Jv(u)  I  i'J,,(Q)- sin(vvj/)  v  =  odd 


(20) 


For  the  transverse  modes  it  can  be  shown  following  the  same  procedure  as  for  the  hybrid 
modes  that 


=  +27ci 

Mj,  =  -27ri 


j,(u) 

J,(u) 


•J,(Q)sin(vi/) 

•J,(Q)cos(v|/) 


(21) 


M,=0 

where  the  values  of  Q  =  r  q  and  v|/  are  the  same  as  defined  before 

Q  =  (a^  +  b^  -  2ab  cos(<j)))’  ^  =  r'-q  and  tg(vj/)  = - kn^i  sin(0)  sin(<j)) - 

Kg  cos(0b)  -  kn^,  sin(0)cos(<j)) 


where 
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q  =  (Kg  008^(83)  +  k^n^,  sin^(0)  -  2kn^,KB  sin(0) 008(615 )cos((t>))’'^ 

The  behavior  of  the  angle  vj/  a8  a  function  of  0  and  <|)  can  be  seen  in  Fig.  3 
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Integral  over  the  variable  r’ 

Now  we  will  calculate  the  third  integral,  i.e.,  we  will  integrate  M  over  the  variable 
r’.  Since  each  component  of  the  vector  M  has  a  different  dependence  on  r’,  each 
component  must  be  calculated  separately.  All  of  the  integrals  are  of  a  similar  type  and  the 
following  relations  will  be  useful 

f  R 

Jr  J.(a  r).J„(b.r).dr  =  ^— ^{b.J„(a  R).J„.,(b.R)-a.J„(b.R).J„_,(a  R)} 

0  a  ~  D 

(23) 

and 

J-Jx)  =  (-!)”  JJx)  (24) 

then,  it  is  obtained  for  the  hybrid  modes 


Mr”  =  27r  •  i 


v-l 


a,  •  R  •  cos(v  -  l)\j/ 


3-'' 


q  •  Jv-,  (u)  V®  -  •  Jv-,  (qR)--^^-^  ' 


Jv(u)  R 


Jv(u) 


+ 


a,  •  R  •  cos(v  +  l)v|/ 


u 


q-Jv.:(u)^^-|-J..,(qR) 

Jv(u)  R 


(25) 


Mr"  =  27ri''-* 


a,  •  R  ■  sin(v  -  l)\j/ 


Jv(u)  R 


Jv(u) 


+ 


+■ 


aj  •  R  •  sin(v  +  l)v|/ 


qJv.,(u)^f^-~J..,(qR) 

Jv(u)  R 


(26) 
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Mp  =  27t  ■  i 


v-l 


a,  •R-sin(v-l)\i/ ^  ^  ^ 


-1 

rJ 


q'Jv-i(u)- 


Jv(u)  R 


Jv(u) 


U  1  2 

-q 


Jv(u)  R 


(27) 


M!“  =27r-i"-' 


a,-R-cos(v-l)y|  u  ( 

r  Jv(u)  R  Jv(u)  J 


q^ 


a,-R-cos(v+l)vi/ j  ^  T  rnP^l 

+ - ^  iq'Jv+iW  ,  X  X  Jv+iiq^)! 


Jv(u)  R 


(28) 


u/B  f  ,  ,  u  ,  ,  Jv.-i(u)l  f-i''  ’cos(v0)  V  =  even 

M,  =  27t— — 2 - jq-J,-,(qR)-— Jv(qR)  ^-  --;  \\ 

2  I  ^  Jv(u)  J  [  1  sin(ve)  V  =  odd 

sJ 


and  in  a  similar  way  for  the  transverse  modes 


(29) 


M.  =.2™-R^  q.J.(qR)-ij,(qR)^ 

R  Ji(u)J 


''uV 


My  =  -27ti 


=  f,.J.(qR)-:^J,(qR)^' 

R  J,(u)J 


(30) 


M,  =  0 

From  the  previous  equations,  it  can  be  seen  that  every  component  of  the  Vector 
Potential  A  has  a  probable  singularity  when  q(0,(|)=O)=u(P)/R.  For  the  hybrid  mode  case 
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v=l  the  previous  condition  gives  the  maximum  amount  of  scatter.  For  other  modes  this 
condition  gives  the  minimum  amount  of  scatter.  Therefore,  for  the  v  =  1  hybrid  mode  or 
for  single  mode  fibers,  q(0,(f)=O)=u(P)/R  gives  the  second  Bragg  condition  for  the 
fabrication  of  a  tapped  Bragg  grating 

kn^,  sin(05j3„ )  =  Kg  cos(0b)  ±  Second  Bragg  Condition  (31) 

Eqs.  16  and  3 1  completely  specify  the  Bragg  parameter  (Kb,  0b)  when  the 
propagation  constant  P  and  the  desired  scattering  angle  0scatt  are  given.  It  can  be  shown 
that  Eq.  3 1  is  a  statement  of  the  conservation  of  momentum  along  a  direction  normal  to 
the  fiber  and  in  the  plane  containing  the  Bragg  vector.  The  ambiguity  in  Eq.  3 1  (the  ± 
ambiguity)  can  be  removed  by  eliminating  the  last  term  or  by  working  in  the  weak 
guidance  approximation  u  ->  0.  Eqs,  16  and  31  reduce  to 


Equations  32  are  the  basic  equations  that  must  be  satisfied  to  fabricate  a  Bragg 
grating.  Fig.  4  shows  what  the  Bragg  period  Ab  should  be  in  order  to  maximize  the 
amount  of  scattered  radiation  in  a  particular  scattering  direction  0scatt  (Oscatt  is  measured 
relative  to  the  fiber  axis  direction).  The  Bragg  angle  0b  is  simply  half  of  the  scatter  angle 
as  shown  in  Eq.  32.  This  result  answers  the  first  question  that  was  addressed  in  the 
introduction,  i.e.,  what  must  the  Bragg  angle  and  period  (0b  and  Ab)  be  in  order  to 


0  30  60  90  120  150  180 


Scatter  Angle  (0  deg) 

Fig.  4:  Bragg  period  as  a  function  scattering  angle  when 
ncl  =  1 .44  and  for  two  different  wavelengths 
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effectively  scatter  radiation  out  of  the  fiber  at  a  desired  scattering  angle  Gscatt-  In  the  weak 
guidance  approximation  ( u  ^  0 ),  the  expressions  for  M  can  be  simplified  by  using  the 
following  relations 

J  (x  0)  =  —  -  ;  K  ](x  ^  0)  =  -  "f"!  i  ~  (3^) 

2  -v!  2  Vx/  x 

in  this  approximation  the  parameters  bi,  b2,  F,  ai  and  a2  converge  to  the  following 
expressions 


V  V  u^  u^  u^ 

bj^— ;  b2  — > - 7;  F->l  +  2— ;  a,— ;  a2->l  +  — 

U  W^  V  V  V 


(34) 


by  replacing  Eqs.  33  and  34  into  Eqs.  25  -  28  the  following  approximations  are  obtained 


cos(v-l)vi/  ^  _qR 

V 


Mf”  =27r  i''  I  ---V  -  ^  -/T  ,  2 .  27tR^vcos(v-  l)x|/ 


=27t-i’ 


sin(v-l)M/  «  ,  qR 


2  •  V— Jy(qR) «  27tR  vsin(v -  l)v|/ 


Jv(qR) 

qR 

Jv(qR) 

qR 


(35) 


=  271  •  i''-’  .  2 .  J^(qR)  »  27TRVsin(v -  1)\|/  • 

q  V  qR 


odd  •  v-1  cos(v-l)\i/  -  ^  qR^  Jv(qR) 


=27t-i 


M  =0 


.2.V— J^(qR)«27tR  vcos(v-l)v|/' 
q  V  qR 


we  have  only  retained  the  terms  of  the  form  v'\  For  the  transverse  modes  it  can  be  easily 
shown  that 


M  =  i  •  27tR^  •  sin(\|/)  ^ 

"  qR 

M  =  i  -  27rR^  •  cos(\|/)-^^^^^ 
"  qR 


M^  =  0 


(36) 
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III.  TBG  BAND  PASS 

For  the  case  of  back  reflection  (  0b  =  90°  and  0scatt=  180°)  it  can  be  shown  (by 
forcing  the  argument  of  the  sine  of  Eq.  1 5  to  be  ±  7i/2)  that  the  spectral  width  of  the 
reflected  energy  is  given  by 


AX  = 


2-N 


(37) 


where  Xo  is  the  wavelength  in  vacuum,  and  N  is  the  number  of  lines  in  the  grating  defined 
from  L  =  N  •  Ab.  Clearly  as  the  number  of  lines  in  the  Bragg  grating  increase,  the  width  of 
the  reflected  component  narrows.  For  a  grating  with  N  =  10,000  lines,  the  reflected  width 
for  a  beam  with  wavelength  of  Xo  =1  pm  would  give  AX  =  0. 1  nm.  Such  a  grating  would 
only  measure  L  =  N  •  Ab  =  N  •  Xo/(2  •  rid)  =  3  mm  (assuming  that  iid  =  1 .5). 

One  of  the  great  advantages  of  Bragg  gratings  is  that,  due  to  the  extreme  spectral 
narrowness  of  the  reflected  component*  (on  the  order  of  AX  =  0  . 1  nm),  many  Bragg 
gratings  can  be  placed  in  a  single  fiber  with  the  possibility  of  interrogating  each  grating 
independantly.  This  can  be  done  by  fabricating  each  grating  in  the  fiber  with  a  different 
Bragg  period  Ab  and  then  interrogating  all  the  gratings  with  a  broad  band  signal  (BBS).  If 
the  spectral  width  of  the  BBS  signal  is  AXbbs  =  50  nm,  centered  at  X  =  1300  nm  in  a  fiber 
of  index  n  =  1.4  and  the  desired  strains  to  measured  are  of  the  order  of  s  =  1,000  pstrains, 
then  the  number  of  sensors  that  could  be  independently  addressed  would  be 
Nsensor=AXBBs/(2nXs)  =  50/3.6  =13  sensors.  Notice  that  in  this  case  the  reflected 
spectral  spread  for  the  desired  strain  measurements  is  AXstrain  =  2nX8  =  3.6  nm  which  is 
much  larger  than  the  spectral  width  of  the  reflected  signal  AX  =  0. 1  nm.  Therefore  it  is 
this  parameter,  AXstrain  ,  that  dictates  the  number  of  sensors  that  will  fit  in  a  fiber  without 
overlap  as  opposed  to  AX .  The  ability  of  placing  many  sensors  in  a  single  optical  fiber  and 
of  being  able  to  interrogate  each  sensor  independently  has  many  potential  applications, 
specially  in  the  area  of  smart  materials  where  the  stress  sta^e  of  an  entire  panel  can  be 
monitored  with  a  single  optical  fiber. 


The  spectral  width  of  the  scattered  energy  when  the  Bragg  angle  is  smaller  than 
0B =1 80°  is  much  larger.  The  basic  reason  for  this  is  that  different  wavelength  can 
effectively  scatter  out  of  a  fiber  at  different  angles,  i.e.,  the  scattering  angle  is  not  fixed  as 
happened  in  the  backscattering  case  described  above.  The  spectral  width  for  this  more 
general  case  will  be  discussed  later  in  this  section.  For  now  we  will  calculate  the  angular 
spread  of  a  monochromatic  (fixed  wavelength  X)  beam  reflected  out  of  a  fiber.  This  can 
be  easily  obtained  by  forcing  the  argument  of  Eq.  15  to  be  equal  to  ±  k/2,  then 


A0 


n^i  •L  sin(0) 


(38) 


14 


NAWCADWAR-95027-4.3 


notice  that  the  beam  can  be  extremely  narrow  as  L,  the  length  of  the  grating,  increases.  In 
general  the  limiting  factor  defining  the  angular  spread  of  the  beam  is  not  the  length  L  of 
the  grating  as  we  have  just  shown,  but  the  spectral  content  AX  of  the  beam,  as  will  be 
shown  next. 


In  the  case  of  a  TBG  the  situation  is  different  from  above  because  now  for  a  set  of 
Bragg  parameters  (Kb  and  0b)  there  will  be  a  broad  range  of  propagation  constants  P  (and 
therefore  wavelengths)  that  will  satisfy  Eq.  16  for  different  scattering  angles.  In  the  weak 
guidance  approximation  ( u->0  or  P  ->  knd )  the  expression  for  the  spectral  width  (AX)  in 
terms  of  the  angular  spread  A0  can  be  obtained  by  differentiation  (dX/d0)  of  Eq.  16 


AX  =  X- 


A0 


(39) 


The  angular  spread,  A0,  can  be  calculated  by  determining  the  angle  at  which  Eq.  35  or 
36  has  its  first  zero  for  (j)=0.  In  the  weak  guidance  approximation  the  function  A(0,(|))  can 
be  written  as 


A(0,<t)) 


cos(0B  )  -  2kn,,  sin(0)) 

q-R 


(40) 


The  zeros  of  the  previous  equation  are  the  zeros  of  the  nth  order  Bessel  function. 
Therefore  we  get 


qR  =  a„ 

sin(0,„„)-sin(0,„„)  =  — ^5 - 

Rkn,, 

where  an  =  2.40,  3.83,  5.14,  6.38  ....  for  n  =  0,1,2,3, ....  are  the  zeros  of  the  Bessel 
function. 

Fig.  5  shows  the  angular  and  spectral  width  for  two  different  fiber  modes  and  two 
different  scattering  angles  as  a  function  of  the  Rk  product.  From  that  graph  it  is  clear  that 
the  larger  the  Rk  product,  the  smaller  the  diffraction  effects,  and  therefore  the  narrower 
the  angular  and  spectral  broadening.  In  contrast,  the  larger  the  Rk  product  is  the  more 
modes  that  the  fiber  will  be  able  to  sustain  and  therefore  more  intermodal  coupling  effects 
in  the  fiber.  For  single  mode  operation  the  Rk  value  has  to  be  smaller  than  the 
2.45/(nci^  -  nco^)^^.  For  a  typical  10  pm  core  fiber  Rk  has  a  value  of  60.  One  can  see  from 
Fig.  5  bottom  that  the  spectral  width  is  AX  =  0.5  pm  for  light  scattered  out  of  the  fiber  at 
90".  This  is  more  then  3  orders  of  magnitude  bigger  than  in  the  case  of  a  back  reflecting 
Bragg  grating.  The  TBG  band  pass  can  be  made  smaller  by  increasing  the  scattering  angle 
as  shown  in  Fig.  5  bottom.  All  these  considerations  will  have  to  be  taken  into  account 
when  fabricating  TBG  sensors. 
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Rk 


Rk 


Fig.  5  Top:  Angular  spread  of  a  TBG  as  a  function  of  the  Rk  value  of  a  fiber  for  two  different 
scattering  angles  and  two  different  modes.  Bottom;  Spectral  Spread  of  a  TBG  as  a  function  of 
the  Rk  value  of  a  fiber  for  two  different  scattering  angles  and  two  different  modes. 
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IV  POYNTING  VECTOR  CALCULATION 


Finally  the  quantity  M  can  be  written  in  the  weak  guidance  approximation  as 


=  a  •  27iR^ v(l  +  v) cos(v -  l)v|/  •  •  5(-Kb  sin(0B )  +  P  -  kn^,  cos(0)) 

qR 

M  (0,(1))  =  a-27iR^v(l  +  v)sin(v- 1)\|/ sin(0B)  +  P- kn  , cos(0)) 

qR 


(42) 


where  “a”  is  a  new  dummy  proportionality  constant  that  contains  various  fiber  parameters. 


From  the  components  of  M  one  can  calculate  the  Poynting  vector  S  and  the  total  radiated 
power 


£lik.(l^y'^|rxA|'f: 

2  Vo  '  ' 


St't-'I 


mV+ImJ'}? 


(43) 


P,«.= 


c%Vn  ,  .8 


1  2  27171  ^  ^ 

x  A|“dn  =  I {|Mer  +|M,f  }sin(e)d8d4.  (44) 

M.  L  3271  e.  J/t'  ) 


where  Me  and  Mj,  are  the  components  of  M  in  a  spherical  coordinate  system.  These 
components  can  be  easily  calculated  fi'om  the  cartesian  components  of  the  vector  M  as 
shown  below 


Mg  =  M,^  cos(0)cos((l))  +  My  cos(0)sin((|)) 
M^  =  -M,j  sin((|))  +  My  cos(<|)) 


(45) 


and  since  Mx  =  M  sin(v-l)ii/  and  My  =  M  cos(v-l)vi/  as  seen  from  Eqs.  42,  then 

M^  +  M^  =  M'  +  My  -  sin"(0)  •  (m,,  cos(<t))  +  My  sin((t)))' 

=  M' (1  -  sin^(0)  sin  ^  ((v  -  l)v|/ +  (!>)) 

therefore 


S(0.‘t>)  =  (1  -  sin^  (0)  sin  ^  ((v  -  \)w  +  ^))-v  (47) 

3271  r  p„ 

In  general,  the  polarization  effects  (which  are  included  in  the  polarization  factor  (1  - 
sin^(0)  sin^((v-l)vi/+<|))))  are  a  much  slowly  varying  fiinction  of  the  scattering  angles  0  and 
(})  than  the  function  M.  As  a  result,  Eqs.  16,  31,  32,  39  and  41  will  remain  valid  even 


17 


NAWCADWAR-95027-4.3 


though  they  were  calculate  directly  from  M  (which  is  proportional  to  the  vector  potential 
A)  rather  than  from  the  Poynting  vector  S  which  has  polarization  information. 

The  total  power  radiated  can  be  calculated  from 


'rad 


3271"  £ 


27171 

J  j  {m'  “  sin"  (0)  sin '  (( V  -  l)vi/  +  <|)))}  sin(0)ded(|)  (48) 


where  in  the  weak  guidence  approximation 


M(0,(t))  =  a  •  27rR"v(l  +  ■  6(-Kj,  sin(0B)  +  P  -  kn^,  cos(0)) 

qR 


(49) 


V.  ENERGY  DISTRIBUTION  OUT  OF  A  TBG 

The  distribution  of  the  energy  out  of  a  TBG  can  be  calculated  from  Eqs.  47  and 
49.  From  those  equations  it  is  possible  to  identify  three  parts  which  contain  the  scattering 
angles  0  and  <)).  The  first  part,  [J(qR)/qR]"  ,  is  a  rapidly  varying  function  of  the  scattering 
angles  and  exhibits  a  singular  behavior  near  qR=0.  This  contribution  is  the  most  important 
contribution  when  determining  the  energy  distribution  out  of  a  TBG.  The  second  part,  the 
delta  function,  places  a  restriction  to  the  values  that  the  polar  scatter  angle  0  can  have  (Eq. 
16).  This  restriction  can  be  directly  incorporated  into  the  first  part.  The  third  part  is  the 
polarization  factor  (1  -  sin"(0)  sin"((v-l)\|/+<|))))  where  vj/  was  defined  in  Eq.  18.  Fig.  6a 
shows  3D  image  of  the  polarization  factor  while  Figs.  6b  and  6c  show  two  cross  sectional 
views  in  parametric  format  of  the  same  polarization  factor  for  a  mode  given  by  v  =  1 .  It 
can  clearly  be  seen  from  these  Figures  that  the  polarization  factor  is  a  slowly  varying 
function  of  the  scattering  angles.  Consequently,  this  contribution  will  not  be  important 
when  determining  the  energy  distribution  out  of  a  TBG. 

Fig.  7a  shows  the  first  part,  [J(qR)/qR]^ ,  of  the  energy  distribution  out  of  a  TBG 
for  two  different  fiber  modes,  the  HEn  and  the  TEoi  and  for  a  scatter  angle  of  0scatt  =  90°. 
Fig  7b  gives  a  pictorial  representation  of  how  the  energy  is  distributed  inside  de  fiber  for 
the  same  two  fiber  modes.  Note  how  the  node  in  the  fiber  for  the  TEoi  mode  translates 
into  a  butterfly  shaped  distribution  out  of  the  fiber.  Other  modes  have  more  complicated 
energy  distributions.  The  HEn  mode  is  the  most  effective  mode  for  coupling  energy  out 
of  an  optical  fiber. 

Fig.  8  shows  the  energy  distribution  for  different  polar  scatter  angles  ranging  from 
0  =  0°  to  0  =  180°  for  (|)  =  0°.  Fig.  9  shows  the  azimuthal  distribution  of  energy  scattered 
out  of  a  TBG  for  various  scatter  angles. 
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Fig.  6a.  Shows  a  3D  parametric  plot  of  polarization  factor  for  the  mode  v  =  1 .  6b  Shows 
a  parametric  plot  of  the  polarization  factor  for  (|)  =  0,  the  vertical  axis  corresponds  to  the 
z-axis,  while  the  horizontal  axis  correspondes  to  the  x-axis.  This  is  the  polar  plane.  6c. 
Shows  a  parametric  plot  of  the  polarization  factor  for  0  =  90°,  the  vertical  axis 
corresponds  to  the  y-axis,  while  the  horizontal  axis  correspondes  to  the  x-axis.  This  is 
the  azimuthal  plane. 
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Fig.  7:  shows  the  first  part,  [J(qR)/qR]^ ,  of  the  energy  distribution  out  of  a  TBG  for  two 
different  fiber  modes,  the  HE, ,  and  the  TEq,  and  for  a  scatter  angle  of  =  90“  Fig. 

7b  gives  a  pictorial  representation  of  how  the  energy  is  distributed  insideThe  fiber  for  the 
same  two  fiber  modes  as  before. 
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Fig.  8.  Energy  distribution  function  for  various  polar  scatterr  angles  ranging  from  0  =  0' 
to  0  =  180° 
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Fig.  9.  Azimuthal  energy  distribution  curves  for  various  scattering  angels. 
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VI.  CONCLUSION 

The  basic  model  and  its  limitations  for  the  fabrication  of  an  optical  fiber  Tapped 
Bragg  Grating  has  been  amply  described  in  section  II  of  this  work.  The  basic  equations 
required  to  fabricate  a  Tapped  Bragg  Grating  (Eqs.  16  and  31)  have  been  derived.  Eqs. 
32,  which  are  the  equations  normally  found  in  the  literature,  are  a  good  approximation  for 
single  mode  fibers  or  for  fiber  that  have  a  small  Rk  value.  The  angular  and  spectral  spread 
of  a  TBG  has  been  derived  in  the  weak  guidance  approximation  (Eqs.  39  and  41).  From 
those  results  it  was  found  that  the  spectral  width  of  a  TBG  is  a  strongly  dependent 
function  of  the  scattering  angle  Gscatt  •  It  was  also  found  that  the  spectral  width  of  a  TBG 
was  several  orders  of  magnitude  larger  than  the  spectral  width  of  a  standard  Bragg 
Grating  back  reflector.  This  results  significantly  limits  the  number  of  TBG  that  can  be 
place  in  a  single  fiber  and  still  be  able  to  independently  interrogate  them.  Finally,  the 
spatial  distribution  of  energy  out  of  a  fiber  has  been  derived  (Eqs.  47).  This  result  is 
important  when  using  TBG  for  demodulation  purposes  or  for  sensor  development 
purposes.  In  particular,  these  results  are  being  applied  for  the  development  of  a  corrosion 
sensor  in  our  group  at  the  NAWCADWAR^. 
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